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Q\ • Electron paramagnetic resonance of ensembles of phosphorus donors in silicon has been detected 

electrically with externally applied magnetic fields lower than 200 G. Because the spin Hamilto- 
nian was dominated by the contact hyperfine term rather than by the Zeeman terms at such low 
£NJ . magnetic fields, superposition states a [ j j) + (3 |||) and — f3 |TI) + a |J. j) were formed between phos- 

phorus electron and nuclear spins, and electron paramagnetic resonance transitions between these 
superposition states and |||) or J.J.) states are observed clearly. A continuous change of a and (3 
with the magnetic field was observed with a behavior fully consistent with theory of phosphorus 
donors in silicon. 
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CN ' PACS numbers: 76.90.+d, 72.20.Jv, 71.55.-i, 76.30.- 
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I. INTRODUCTION 

A phosphorus in silicon is attracting much attention 
towards realization of solid-state quantum information 
processors. It can be viewed as a two-qubit system 
having one 31 P nuclear spin (/ = 1/2) and one elec- 
tron spin (S = 1/2) J^^^ Coherent manipulation of 
its electronic states^^ 9 - coherent transfer of states be- 
tween electron and nuclear-spins,— and large hyperpo- 
larization of nuclear spin s 11 ' 12 have been demonstrated 
recently. Thanks to phosphorus's long spin dephasing 
time ) 13 ' 14 enrichment of silicon with nuclear spin-free 
28 Si has suppressed the background isotope fluctuation 
significantl y 15 ' 16 to make possible the optical detection 
of 31 P nuclear spin states*^ Because a standard elec- 
tron paramagnetic resonance (EPR) measurement re- 
quires at least 10 9 spins or more, much more sensitive 
electrical detection methods of phosphorus EPR have 
been attracting attention^^O^^^^ 24 * 2 ^ The 
record so far reported is the detection of ~ 50 phos- 
phorus spin states 2 ^ and extensive efforts are underway 
worldwide to detect single phosphorus spin states^ All 
of the previous phosphorus EPR studies have been per- 
formed in the "high-magnetic-field regime," which can 
be defined as B 3> 200 G for phosphorus in silicon as 
we demonstrate later, and have observed two EPR al- 
lowed transitions i 6 -^^ ' 11 ' 19 ' 24 ' 25 ' 26 The present work re- 
ports electrically detected magnetic resonance (EDMR) 
of phosphorus spin states in silicon which shows five of six 
possible transitions expected for the phosphorus donors 
in silicon under the low magnetic field B < 200 G. 

The spin Hamiltonian of an isolated phosphorus atom 
placed in an externally applied magnetic field B is given 
by; 



Hsup = 9eHBBS z - g n [i n Bl z +aS ■ I, 



(1) 



resent the electron Zeeman, nuclear Zeeman, and con- 
tact hyperfine interaction between phosphorus electron 
and nuclear spins, respectively. Here g e iiB/2nh ~ 28 
GHz/T and g n ii n /2ith w 17.2 MHz/T are given by elec- 
tron and nuclear ^-factors g e sw 1.9985 and g n « 2.2632, 
respectively 28 ' 29 The hyperfine constant is a/2Trh as 
117.5 MHz. 3 - Eigenstates of this spin Hamiltonian are 
given by; 



|i> = ITT), 

|2) = a\U)+P\m, 

|3) = -/3|U)+«|IT), 

|4) = |U>, 



(2) 
(3) 
(4) 
(5) 



where a — cos | and (3 — sin ^ . 77 is the angle be- 
tween externally applied magnetic field direction and ac- 
tual electron and nuclear spins precession axis given by 
tan 77 = A magnetic quantum number 



where S and J are electron and phosphorus nuclear spins, 
respectively. The first, second, and third terms rep- 



+i (-5) is represented by f (!) and an arrow on the left 
(right) in each ket represents the electron (nuclear) spin 
state. Fig. 1(a) shows the magnetic field dependence of 
the four states expected for phosphorus in silicon. At the 
high magnetic fields (B ^> 200 G) a ~ 1 and (3 ~ 0, i.e., 
the four states simply become ||T), ITI), I IT), an( l III)- 
The EPR allowed transitions in this regime are limited to 
two: ITT) «■ III) and |t|) |||). In the low-magnetic- 
field regime defined by B < 200 G, the eigenstates of 
1 2) and |3) change continuously because r], i.e., a and (3, 
change significantly with B as shown in Fig. 1 (b) . There- 
fore, the degree of superposition between ||J.) and |||) 
that determines the EPR allowed transitions also changes 
with B. For example, the transition |1) |2) at the high 
magnetic field corresponds to the nuclear magnetic res- 
onance (NMR) HI) and it cannot be observed 
as the EPR. However, the same transition becomes EPR 
observable at the low magnetic field because the EPR al- 
lowed component ||T) & III) emerges with (3. Note that 
transitions |2) |3) and |1) <^> |4) are weak because 
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FIG. 1: (a)Externally applied magnetic field dependence of 
the spin states of phosphorus in silicon defined by Eqs. ([2} - 
©. Six allowed transitions |1) O |3), |2) O |4)C |2) ^ |3>, 
1) O |4), |1) <^ |2), and |3) |4) are labeled by □, ■ 
A, A. •, and o, respectively, (b) Externally applied magnetic 
field dependence of a and (3. a —> 1 and f3 — > for B 3> 200 
G. 



they are allowed only in the second order. Nevertheless, 
two aspects of quantum control that cannot be realized 
in the high-magnetic-field regime is expected to become 
possible in the low-magnetic-field regime; 1) controlling 
the ratio of a and (3 to change the degree of superpo- 
sition by the magnetic field and 2) changing the popu- 
lation of the four states by utilizing the six transitions 
that are made allowed. The present work demonstrates 
these properties experimentally using electrical detection 
of phosphorus EPR and develops a quantitative theoret- 
ical model to support our observation. 



II. EXPERIMENTAL 

A sample was a bulk Czochralski-grown n-type silicon 
single crystal having phosphorus concentration ~ 10 16 
cm~ 3 . It was cut into a rectangular shape of the di- 
mension 8x2x1 mm 3 . Ohmic contacts were prepared 
at both ends of the long axis by arsenic implantation of 
2 x 10 15 cm~ 2 at 25 keV followed by annealing at 980 °C 
for 25 seconds and vacuum deposition of the 5-nm-thick 
palladium and 50-nm-thick gold layers. The sample was 
placed in a cryostat with optical windows. A white light 
from a halogen lamp placed outside of the cryostat was 
focused onto the sample through the optical window for 
steady state excitation of the electron-hole pairs to main- 
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FIG. 2: EDMR signals (change in the sample photoconductiv- 
ity under continuous white light illumination from a halogen 
lamp) vs. externally applied magnetic field under irradiation 
of different 500 mW radio frequencies (RF) as indicated in 
the figure. The sample is phosphorus-doped bulk silicon sin- 
gle crystal ([P]~ 10 16 cm" 3 ) kept at T = 5 K during the 
measurement. The peaks indicated by □ correspond to the 
transition jl) <4> |3) labeled by the same mark in Fig. []Ja). 
Likewise, ■ A, •, and o correspond to |2) <4> |4)C |2) <4> |3), 
|1) <4> |2), and |3) O |4), respectively. Peaks labeled by * 
correspond predominantly to the paramagnetic resonance of 
the interface center for the reason discussed in the text. 



tain the sample resistance at ^10 kJl . The sample was 
connected with a series resistor of 10 ktt. A constant 
voltage of typically 10 V was applied to the series of the 
sample and resistor. A coaxial cable was used to connect 
a RF source with an irradiation coil whose opposite side 
was connected to a 50 terminator. Externally applied 
magnetic field was provided by a 300 mm bore electri- 
cal magnet. Another pair coil was placed in the cryostat 
to modulate the externally applied magnetic field for the 
lock-in detection of the divider voltage corresponding to 
the change in the sample photoconductivity (EDMR sig- 
nal). 



III. RESULTS AND DISCUSSIONS 
A. EDMR peak positions 

Fig. [2] shows the EDMR signals obtained at T = 5 
K. By changing the irradiation frequencies, six different 
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FIG. 3: A plot of irradiated frequency vs. externally ap- 
plied magnetic field showing positions of experimentally de- 
termined peaks represented by the same marks as in Figs. [1] 
and The marks indicated by □, ■. A, •, o, and * cor- 
respond to |1) o |3), |2) |4)C |2) o |3), |1) |2), 
3) j4), and interface center transitions, respectively. Solid 
curves are rigorous theoretical calculations of resonance posi- 
tion of phosphorus in silicon, which show excellent agreement 
with experimentsD 



transition peaks labeled □, ■. A, •, o, and * are observed 
clearly. 

Fig. [3J shows externally applied magnetic field vs. 
RF frequencies of the six observed resonance positions. 
Solid curves are theoretically expected results for differ- 
ent transitions of phosphorus in silicon using Eqs. ([2]) - 
(O that are calculated rigorously with no fitting param- 
eter. The excellent quantitative agreement between our 
experiment and theory supports the successful observa- 
tion of the transitions listed in the captions of Figs. [TJ [21 
and [3] 

The transition indicated by * agrees very well with the 
theoretically predicted resonance positions for |1) 4=> |4) 
labeled by ▲. However, we conclude that this peak is 
composed predominantly by the EPR transition of the 
paramagnetic defects situating around the interface be- 
tween bulk Si and native surface oxide SiC>2 for the follow- 
ing reasons. As apparent from Fig. [TJa), the transition 
|1) |4) indicated by ▲ corresponds to a typical elec- 
tron Zeeman transition whose energy is approximately 
proportional to B. Therefore, EPR transition energies of 
other defects with very little hyperfine interaction, such 
as those of the interface spin states ; 32 ' 33 overlap with 
this transition. Here the EPR of the interface spin state 
is allowed in the first order and, therefore, expected to 
have strong intensity. However, the transition |1) |4) 
labeled by A is allowed only in the second order and 
should demonstrate approximately the same intensity as 
the transition |2) <^ |3) labeled by A in Fig. 2. The fact 
that the intensity of * being much stronger than that of A 
suggests that the * transition arises predominantly from 
the EPR transition of the interface defects that is allowed 
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FIG. 4: (a) A schematic diagram of the EDMR mechanism. 
An electron at the phosphorus level undergoes spin resonance 
and falls to the interface level when the spin direction of phos- 
phorus state and that of interface state form a spin singlet. 
Once the electron bound to phosphorus is gone, the phos- 
phorus capturers another electron from the conduction band, 
leading to the change in the photoconductivity. Here the con- 
duction electrons are captured by phosphorus at the rate G 
and electrons at phosphorus go back to the conduction band 
at the rate D or are captured by the interface states at the 
rate R. (b) RF power dependence of the EDMR signal for 
the transition |2) <4> |4). The solid curve is the fitting using 
Eqs. rtT2jl and (ITS"!) . 



in the first order. The importance of the presence of the 
interface defects was confirmed by removing the surface 
oxide using a dilute hydrofluoric (HF) solution. The HF 
treatment made all of the EDMR signals nearly invisi- 
ble. Leaving the sample for a few days in air to cover 
the sample surface with native oxide again retrieved the 
intensity of the all EDMR signals completely. 



B. RF power and magnetic field dependencies of 
the phosphorus EDMR 

The fact that the presence of the surface oxide is 
needed to observe the EDMR signals allows us to develop 
a theoretical model describing the RF power and mag- 
netic field dependencies of the EDMR signal intensity 
as the following. We assume that the interface defects 
act as spin-dependent recombination centers for electrons 
bound to phosphorus donors, fn this sense, our model 
is an extension of existing two-level spin-dependent re- 
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TABLE I: Recombination rate of each transition between 
phosphorus and interface spin states.— 



transition. In this case; 
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combination j^de^SOMSULSi^MLmm^ to the 
four levels. The spin Hamiltonian of the system during 
EDMR is given by 

Hedmr = g e ^BB S z + g e ^ B BiS x cos (ut) + aS ■ I. (6) 

Here the nuclear Zeeman term is neglected because g n ^, n 
is approximately 10 3 times smaller than g e fJ-B and the J 
coupling term (JS ■ Si where S and Si are the electron 
spins of the phosphorus and interface states, respectively) 
is also neglected assuming a^> J. The second term arises 
from the RF irradiation and this perturbation term is 
defined as 7i to calculate the transition probability W 
for the system described by Eq. ^ using the Fermi's 
golden rule, 



W = ^\{f\H'\i)\ 2 6(E f -E i -hv). (7) 



The electron-hole recombination via interface states takes 
place only when S and Si form a spin singlet and does 
not occur when they form a spin triplet to establish a 
"spin blockade."— 

Other important essences of our model are described 
in the caption of Fig. 0Ja). Let us consider \i) where 
i = 1, 2, 3, or 4 is one of the four phosphorus spin states 
as defined by Eqs. (|2l) - © and \a) where a = | or 
I corresponds to spin up or down of the interface state, 
respectively. Using G, D, and R defined in Fig. HJa), we 
obtain the rate equation 



d 

dt Ni > 



G(N - 



jo-' 



(D + R ia )N io 



(8) 



where N is the total number of electron pairs and Ni a is 
the number of electron pairs in spin states i and a. The 
recombination rates Ri a that have been obtained using 
the method described in Ref. 43 are listed in Table fl] 
The pairs with (i, a) with Ri a = correspond to the spin 
blockade. Now we consider a representing example where 
the irradiated RF is in resonance with the |2) <^> |4) 



d_ 

dt 



N 2a = G\ N-J2 N ^' 

3<r' 

W (N 2a - N ia ) , 



(D + R 2a ) N 2a 



= G IN - J2 N 3°> - ( D + R^) N io 



+W{N 2a -N ia ), 



(9) 



(10) 



where W oc (g e ^BBi/2) 2 a 2 is the transition probability 
that is proportional to the RF irradiation power around 
the origin4i This partly lifts the spin blockade and pro- 
motes the recombination. The steady state solution for 
the number of electron recombining I(W) becomes; 



I(W) 



NG 



! + <?[£ 



i(T D+R ia 
Ri< 



E 



X(W)=J2 



D + R t 
1 



WX(W) 
- + WDX(W) 



(11) 



(D + R 2a ) (D + R i0 

(R2a — Ria) 2 



(12) 



(2D + R 2a + R 4 a) W+(D + R 2a ) (D + R 4a ) 
Then the EDMR signal intensity Sedmr is obtained as; 
Sedmr = I(W) - 1(0) 

= NG D + 8G 



1 + G E«T D + Ria 

WX(W) 



1 + G 



y. 



D+Ri, 



WX(W) 



(13) 



This result shows that the signal intensity is propor- 
tional to the irradiated RF power and, therefore, to W 
around the origin. This corresponds to our experimental 
observation shown in Fig. IDJb) when Aa/a is defined as 
Sedmr/ 1(0)- A solid curve shown in the figure is the 
successful fitting by Eqs. ^T2j) and using N, G, D, 
and R as fitting parametersD A set of appropriate values 
we found are N = 6.4xl0 5 cm" 3 , G = 5.0xlCT 9 sec" 1 , 
D = 2.3X10 1 sec" 1 , and R = UxlO 1 sec" 1 . 

Finally, we show in Fig. [5] the magnetic field depen- 
dence of the EDMR signal intensity. We used the RF 
power of 500 mW, which was large enough to saturate 
the signal as shown in Fig. 0Jb). Here If Now 
the two WX(W)'s in Eq. (fL2|) are replaced by a constant; 



WX(W) -> 



1 



(R2a — RiaY 



(D + R 2a ) (D + R ia ) 2D + R 2a + R ia ' 

(14) 

Such relations with appropriate Ri a have been used to 
fit representative experimental results shown in Fig. [S] 
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External Magnetic Field (G) 

FIG. 5: EDMR intensity vs. externally applied magnetic 
field. Experimentally determined positions are represented 
by the same marks in as Fig. [5] The peaks indicated by □, 
■ , A, •, o, and * correspond to |1) <S> |3>, |2) |4)C |2) 
|3), |1) <4> |2), |3) <4> |4), and interface center transitions, 
respectively. Solid curves are fittings for allowed transitions 
□ , ■. •, and o using Eqs. CSJ and JTUl with WX(W) as 
described in the text. The fits for • and o overlap completely. 
The interface center transition * is fitted with Eqs. (|13[) and 

(TSJ. 



This result has been used to fit the behavior of the inter- 
face center peak (*) in Fig. [5] A set of appropriate values 
we found are N = 8xl0 3 cm' 3 , G = 2.0xlCT 5 sec" 1 , D 
= 3.6X10 1 sec' 1 , and R = 3.9X10' 1 sec' 1 . Here the 
values of AT, D, and R are approximately the same as 
the ones obtained for phosphorus but the value of G is 
different. 

A significance of the results shown in Fig. \5\ is that 
the intensity of transitions changes below 200 G in ac- 
cordance with theory. The standard EPR allows for ob- 
servation of only |1) <J4> |3) and |2) <S4> |4) because (3 = 
in the high magnetic fields. However, the value of (3 in- 
creases with decreasing the magnetic field, especially for 
the magnetic fields below 200 G, and approaches 1/V2 
as B — > 0. Naturally, the intensity of 1 1) |3) and 
1 2) 4=> 1 4) decreases because the components of the EPR 
allowed | TT> lit) and ||J.) ^ |U) diminish. For the 
same reason transitions such as |1) <=> |2) and |3) |4) 
appear only when B < 200 G. This observation leads us 
to conclude that it is possible to form the superposition 
states between electron and nuclear spins of phosphorus 
in the regime of the low magnetic field B < 200 G. Their 
superposition coefficients a and (3 can be controlled sim- 
ply by selecting an appropriate magnetic field. 



Note that our model is not applicable to second-order al- 
lowed transitions such as |2) |3) and |1) |4). From 
the fitting of the first-order allowed |2) <^> |4) transi- 
tion, we obtain a set of appropriate values N — l.OxlO 4 
cm" 3 , G = 6.7xl0~ 6 sec" 1 , D = 2.9X10 1 sec' 1 , and R 
— 7.0xl0 _1 sec -1 . Note that the magnetic field depen- 
dence of a and (3 employed here is same as the one shown 
in Fig. [TJb). Therefore, the excellent agreement between 
our experiment and model shows that the coefficients of 
superposition can be controlled by the choice of B and 
reach the maximally entangled states 1/V2 (|T4) + |4T)) 
in the limit of B = as expected4£ 

Similarly, we can derive a relation for the interface cen- 
ter EDMR, since the rotation of the interface electron 
spin Si also lifts the spin blockade and enhances the re- 
combination. We consider again a rate equation using N, 
G, D, and R to obtain the same Sedmr as Eq. (fUSj) but 
different X(W) from the phosphorus resonance case4£ 



1 



^ (D + R^) (D + R tl ) 



{2D + i?. lT + Rn) W + (D + R n ) (D + Rii) 



(15) 



When the signal intensity is saturated, T / FX(PF) is given 
by; 



WX{W) -> 



1 



f(D + Rit)(D + R41) 2D + R ir +R il ' 

(16) 



IV. CONCLUSION 

Electron paramagnetic resonance spectroscopy of an 
ensemble of phosphorus donors in silicon has been per- 
formed successfully to map out the behavior of phos- 
phorus spin states at the magnetic field lower than 200 
G. Formation of the superposition states a |f J.) + [3 |J.f) 
and — /3|T4) + a lit) has been confirmed with the val- 
ues of a and (3 changing continuously with the magnetic 
field in accordance with theory of phosphorus in silicon. 
Dependencies of the EDMR signal intensity on the RF 
power and magnetic field have been described success- 
fully by a model assuming a spin-dependent recombina- 
tion of phosphorus electrons via defects situating around 
the oxide/silicon interface. 
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